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ABSTRACT 

The  modified  phase-space  theory  of  reaction  rates  is  applied  to  the 
dissociation  of  HC&  in  a  heat  bath  of  argon  atoms.  Excellent  agreement 
is  obtained  between  the  theoretical  predictions  and  the  shock-tube  measure¬ 
ments  of  the  dissociation  rate  coefficient  over  the  temperature  range 
2500-5000  °K.  The  recrossing  correction  factor  and  nonequilibrium  correc¬ 
tion  factor  are  obtained  from  Monte  Carlo  trajectory  calculations  for 
states  near  the  dissociation  limit.  The  trajectories  were  sampled  within 
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and  was  used  to  solve  the  steady-state  master  equation.  Important 
features  of  collisions  of  highly  asymmetric  diatomic  molecules  are 
discussed,  and  several  typical  reactive  trajectories  are  shown  to 
Illustrate  the  importance  of  rotational  motion  in  such  collisions. 

1.  INTRODUCTION 

The  modified  phase-space  theory  of  reaction  rates  has  been 
developed  by  Keck  and  his  co-workers*”^  and  applied  quite  success¬ 
fully  to  the  three-body  recombination  and  dissociation  of  diatomic 
molecules  which  have  similar  atomic  masses.  However,  as  demonstrated 
by  Shul,  Appleton,  and  Keck,  the  theoretical  predictions  of  the 
dissociation  rate  coefficients  k^  for  HF  and  HC&  are  much  too  low  by 

comparison  with  experimental  measurements  (see  Figures  9  and  10  of 

B 

Reference  A),  although  the  barrier  rate  coefficients  k,  lie  above  the 

a 

experimental  measurements  and  have  about  the  right  temperature  depen¬ 
dence.  The  reason  for  this  disagreement  was  believed  to  be  that  the 

estimates  of  the  recrosslng  correction  factors  (N/N  )  £  0.01,  and  the 

o 

nonequilibrium  correction  factors  (k/kft)  =  0,01  were  too  small. 

The  correlation  formulae  used  in  Reference  A  to  estimate  the 

values  of  (N/N  )  and  (k/k  )  were  deduced  using  the  results  of  Monte 
o  e 

Carlo  trajectory  calculations^  and  master  equation  solutions^  In  which 

the  recombining  (dissociating)  molecules  had  atoms  of  similar  masses. 

4 

It  was  suggested  that  the  mechanics  of  reacting  collision  processes  in 
which  one  of  the  recombining  atoms  had  a  very  small  mass,  such  as  the 
hydrogen  halides,  was  sufficiently  different  from  those  cases  considered 
in  References  6  and  7,  that  the  correlation  formulae  were  no  longer  valid. 
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Thls  paper  presents  results  of  Monte  Carlo  trajectory  calcula¬ 
tions  for  the  HC&  +  Ar  system,  and  from  these,  determinations  were 
made  of  the  dissociation  rate  coefficients  which  could  be  compared 
with  shock-tube  measurements.  In  Section  II  the  equations  and  statisti¬ 
cal  methods  used  to  obtain  the  numerical  results  are  derived  and  summa¬ 
rized.  Section  III  presents  a  summary  of  the  results  and  their  com¬ 
parison  with  the  experimental  measurements,  and  Section  IV  contains 
our  concluding  discussion. 

II.  MONTE  CARLO  TRAJECTORY  CALCULATIONS 
The  steady  state  reaction  rate  coefficient  k  (subscript  r  for 
recombination;  d  for  dissociation)  is  given  by  the  relationship 

k  -  kB  (N/N  )  (k/k  )  (2.1) 

o  e 

where  the  "barrier"  rate  coefficient  k  is  defined  as 

kB  -  rtv(0)/V  KfMjl  .  (2.2) 

In  Equation  (2.2),  V  is  the  normalization  volume  in  configuration  space, 
[Mj]  is  the  concentration  of  particles  of  type  in  the  Initial  state, 

#y(e)  -  P(v  •  n)  ds  (2.3) 

is  the  "variational  rate"  (see  Reference  2)  for  the  surface  S(c)  defined 
by 

(H  -  B)/kT  -  C  ,  (2.4) 

and  the  corresponding  equilibrium  one-way  crossing  rate&(E)  is  given 
by  the  expression 

^g(c)  -  (N/Nq)  -ftv(e) 

Here,  p  is  the  density  of  the  representative  points,  v  is  the  generalized 
velocity,  n  is  the  unit  outward  normal  to  ds,  and  B  and  H  are  the  barrier 
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helght  and  Hamiltonian  of  the  molecule,  respectively.  As  Implied  by 
Equation  (2.4),  the  energy  (e)  of  the  molecule  is  measured  from  the 
top  of  the  rotational  barrier.  Other  definitions  and  more  detailed 
descriptions  can  be  found  In  References  2  and  3. 

Monte  Carlo  trajectory  calculations  as  applied  here  efficiently 
evaluate  the  correction  factors  (N/Nq)  and  (k/kg)  by  statistically 
sampling  within  the  reaction  zone  and  numerically  Integrating  the 
classical  equations  of  motion  In  both  time-vise  directions  to  deter¬ 
mine  the  complete  history  of  the  collisions.  By  sampling  with  a 
weight  proportional  to  the  local  flow  rate  p(v  *  n)  across  the  surface 
S(e),  a  statistical  distribution  of  trajectories  is  obtained  which 
reflects  their  a  priori  contribution  to  the  reaction  rate,  thus  enabling 
calculations  to  be  made  of  the  correction  factors  without  further  weight¬ 
ing.  The  representative  points  were  selected  from  the  reaction  zone 
region,  l.e.,  on  the  surface  S(e),  and  the  equations  of  motion  were 
Integrated  outwards  In  both  time-wise  directions  and  terminated  as 
soon  as  the  Interaction  with  the  third  body  became  negligible.  This 
avoided  unnecessary  Integrations  In  regions  where  interactions  were 
very  weal  and  inconsequential  to  the  outcome  of  the  collisions. 

IIA.  Recrossing  Correction  Factor 

The  Nq  trajectories  sampled  which  cross  the  barrier  surface, 
e  ■  0,  can  be  divided  Into  four  categories,  l.e.,  f  ■*  b,  f ■  +  f ,  b  ■*  b, 
and  b  f ,  depending  on  whether  the  Initial  and  final  states  of  the 
molecule  are  free  (f)  or  bound  (b).  The  f  ■*  b  trajectories  are 
"reactive"  (In  terms  of  recombination)  while  the  others  are  all  regarded 

as  being  "non- reactive. "  The  recrossing  correction  factor  (N/N  )  Is 

o 
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then  given  by  the  expression^ 

N/Nq  -  Z[N(f|n|b)/n]/No  (2.6) 

n 

where  N  is  the  effective  number  of  reactive  trajectories  and  N(f|n|b) 
is  the  number  of  trajectories  that  go  from  a  free  state  f  to  a  bound 
state  b  with  n  crossings  of  the  barrier  surface  in  the  direction  f  b, 
and  n  -  1  crossings  in  the  direction  b  -*■  f. 

IIB.  Transition  Rates  and  Nonequilibrium  Correction  Factors 
The  nonequilibrium  correction  factor  (k/kfi)  is  the  ratio  between 
the  steady  state  reaction  rate  and  the  corresponding  equilibrium  reac¬ 
tion  rate  coefficient.  Keck  and  Carrier^  Investigated  the  coupled 
vibration-dissociation-recombination  process  for  gas  mixtures  which  con¬ 
sisted  primarily  of  diatomic  molecules  highly  diluted  in  a  heat  bath  of 
inert  gas  atoms  such  as  argon.  They  presented  techniques  for  solving 
the  appropriate  master  equations  and  showed  that  great  simplifications 
could  be  obtained  if  the  kernel  RCe^,  ef)  in  the  master  equation  was 
separable.  The  kernel  R(eit  ef)  is  the  "one-way"  equilibrium  transition 
rate  between  energy  states  e  ■  and  e  »  ef,  per  unit  volume,  per 
unit  initial  and  final  energies.  The  energies  are  measured  in  the  units 
of  kT. 

Monte  Carlo  trajectory  calculations  were  first  used  to  generate  a 
relative  numerical  distribution  of  R(e^,  e^)  as  a  function  of  and  e^. 

A  functional  form  for  Rte^,  ef)  was  then  obtained  by  a  suitable  fitting 
procedure  which  approximated  the  numerical  distribution.  This  functional 
form  was  used,  following  Keck  and  Carrier's  analysis,  ^  to  solve  the 
appropriate  master  equation  and  thus  to  obtain  the  corresponding  non¬ 
equilibrium  correction  factors  k/k^. 
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The  method  of  generating  the  relative  numerical  distributions  of 
RCc^i  )  from  Monte  Carlo  trajectory  calculations  employed  was  similar 

g 

to  that  used  by  Mansbach  and  Keck  which  they  applied  to  the  problem 
of  atomic  excitation  and  ionization  by  thermal  electrons.  Briefly, 
samples  were  taken  from  surfaces  within  a  few  kT  below  the  rotational 
barrier  (e.g. ,  e  ■  -1,  -2,  -3)  in  addition  to  the  samples  taken  on  the 
barrier  surface  (e  ■  0).  The  reactive  trajectories,  i.e.,  f  -*■  b,  were 
divided  into  groups  or  "boxes,"  according  to  the  magnitude  of  and 
,  and  the  effective  number  of  reactive  trajectories  in  each  box  was 
then  counted.  The  equilibrium  transition  rate  is  given  by 


R(eit  ef)  -  [fk(e)/(Aen  -rtv<e) 


(2.7) 


where  (Ac)  is  the  size  of  the  box,  e  is  the  position  of  the  surface 

from  which  samples  were  taken,  the  subscript  k  specifies  the  box  defined 

2 

by  Ei ,  e,  and  (Ae)  ,  and 

N 

o 


f.  (e)  -  Z  (x,. /n,)/N  (e) 
*  j„l  JR  J  o 


(2.8) 


is  the  effective  number  of  reactive  trajectories  falling  in  the  k-th 
box.  Here,  x^k  -  1  if  the  trajectory  is  reactive,  x^k  -  0  if  it  is 
nonreactlve,  and  n^  is  the  number  of  crossings  of  that  trajectory  in 
the  reactive  direction  (e.g.,  from  free  to  bound,  in  the  case  of  recom> 
bination) . 

As  was  observed  previously,  a  separable  form  of  the  kernel  R(e1,  e^) 
results  in  great  simplification  in  solving  the  master  equations.  He 


found  that  the  empirical  form 

R(ei,  ef)  «AG(ef)  exp(o  ef  -  8^) 


(2.9) 


was  able  to  reproduce  the  relative  numerical  distribution  of  R(e^,  e  .) 
with  sufficient  accuracy.  In  Equation  (2.9)  A,  a,  and  @  are  constants 
and  G(Cf)  is  a  relatively  slowly  varylrg  function. 

It  was  observed  from  the  map  of  th»  relative  numerical  distirbu- 
tlon  of  R(e^,  Ef)  that  the  kernel  exhibits  a  sharp  maximum  along  the 
diagonal  ,  dies  off  very  rapidly  away  from  the  diagonal,  and 

varies  relatively  slowly  along  the  diagonal.  Thus  the  fitting  was 
made  easier  by  rotating  the  energy  coordinates  and  expressing  the 
kernel  in  the  form  R(A,  e),  where 


defined  as 


-8- 


Km(e)  "  £>  def  /"de1[R(e1,ef)^(e)]  (eA  -  e^®"1  .  (2.14) 


By  substituting  Equations  (2. 7) and  (2.8)  into  Equation  (2.14)  and 

evaluating  the  integral  by  Monte  Carlo  methods,  we  obtain 
N 


K  (e)  -  Z  (x,/n.)  (e.  -  e£)"“1/N„(e) 
®  j«l  3  j  1  1  J  ° 


(2.15) 


where  ■  1  if  the  trajectory  is  reactive,  and  x^  ■  0  otherwise.  This 
is  the  equation  used  in  the  actual  calculation  of  the  moments  K  (e)  from 
data  generated  by  trajectory  calculations. 

Now,  the  equilibrium  one-way  crossing  rate  ^  (c)  is  defined  as 


ft  <e>  E  JL  def  dei  R(ei»  ef)  *  (2.16) 

Rotating  the  energy  coordinates  to  [A,  e]  and  using  Equation  (2.12),  we 
obtain 

(R,  (e)  sAG(c)  exp(-be)/(a2  -  b2/4)  (2.17) 

where  G(e)  has  been  expanded  about  e  and  small  terms  dropped. 

In  addition,  we  found  by  plotting  (N/Nq)  versus  C,  that  the  form 

N(e)/No(e)  -  C  exp(pe)  (2.18) 

with  both  C  and  p  constant,  correlates  the  numerical  data  reasonably  well. 
By  combining  Equations  (2.5),  (2,17),  and  (2.18),  we  obtain 


(ftv(e)  expCS)/tyO)  -  G(e) 

(2.19) 

b  -  1  -  p  , 

(2.20) 

and 

A/W^O)  -  C(a2  -  b2/4). 

(2.21) 

The  moments  defined  by  Equation  (2.14)  can  then  be  evaluated  by  using  the 
above  relationships  as  follows: 
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K  (e)  -  A[G(e)  expO-eJfctt)]"1  I^J.l  de  /“  dA  G(e)  exp(-aA-be)  A®"1.  (2, 

ID  V  £— il/4  O 

When  a  »  b/2 ,  a  good  approximation  is 


K  (e)  js  C  m!  a^"m  exp[(l-b)e] 
m 


(2.23) 


The  ratios  of  these  moments  are  independent  of  the  energy  e,  and  the 
parameter  a  may  be  obtained  directly  from  any  of  these  ratios.  However, 
since  we  anticipate  from  the  diffusion  analysis^  that  it  is  the  first 
and  the  second  moments  which  are  most  important  for  determining  the 
steady-state  rate  coefficients,  we  have  used  the  value  of  a  given  by 
the  ratio  2^/^. 

Having  obtained  a  separable  form  for  the  transition  kernel,  we  can 
now  proceed  to  solve  the  steady-state  master  equations  which  will  yield 
the  nonequilibrium  correction  factor.  Keck  and  Carrier ^  showed  that 
the  steady-state  master  equation  can  be  transformed  exactly  into  an 
equivalent  diffusion  equation  and  then  solved  subject  to  the  appropriate 
boundary  conditions.  For  the  case  of  the  separable  kernel,  they  assumed 

r«(e.)  r.  (e.)  ;  e  >  e. 


ef) 


(  2V  i/  lv  f'  ’if 
r2(ef)  ri^ei^  ;  ef  >  ei 


(2.24) 


and 


R(c,  e)  -  r^e)  r( 


(2.25) 


where  r  is  a  constant  and 
c 


R(c,  e)  -  /  R(e',  e)  de' 


(2.26) 


is  the  equilibrium  recombination  (or  dissociation)  rate. 


By  solving  the  equivalent  diffusion  equation,  they  obtained  the 
expression  for  the  nonequi librium  correction  factor 

k/ke  -  tffcO)  [f°_s  WZ"2  de  +  r2(0)/(rcZ(0))J}'1  (2.27) 

where<£(0)  is  the  one-way  equilibria  flux  of  molecules  across  the  energy 
surface  e  ■  0,  defined  in  Equation  (2.16), 

W(e)  -  r2(dri/de)  -  r^dr^de)  (2.28) 

is  the  Wronskian  of  r2  and  r^  , 

Z(e)  -  r2(e)  r^e')  de’  +  r^eH/^  r2(e')  de'  +  rcJ  (2.29) 

is  the  equilibrium  collision  rate,  and  6  -  D/kX  is  the  dimensionless 
potential-well  depth  of  the  molecule. 

Since  we  sampled  trajectories  for  recombination,  the  condition, 

Ci  >  cf*  was  satisfied  in  all  cases.  Thus,  for  the  kernel  given  by 
Equation  (2.9),  we  have 

r2*ei^  "  exP(”8ei)  ,  r^Cj)  ■  A  G(ef)  exp(ac£)  (2.30) 

and 

r  -  1/B  .  (2.31) 

c 

By  substituting  Equations  (2.34)  and  (2.35)  into  Equations  (2.32)  and 
(2.33)  and  then  into  Equation  (2.31),  we  obtain 

«  [1  -  (<*/3)2]  [G(e*)/G(0)J 

e 

* 

where  £  is  the  location  of  the  maximum  of  G(e)  exp (-be). 


(2.32) 
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III.  RESULTS 


IIIA.  The  Interaction  Potential 

A  Morse  potential  function  is  used  to  describe  the  two-body  interac- 

4 

tion  between  H  and  Cl  as  before,  with  the  same  potential  parameters; 

o  o 

i.e.,  D  ■  4.62  eV(53600  °K),  rg  ■  1.275  A,  and  3  ■  1.87  A-  .  However, 

3  6 

dumb-bell  models  as  used  before  '  for  the  three-body  interaction  poten¬ 
tial  are  not  well  defined  when  VTO  and  are  considerably  different 
and  have  attractive  potential  wells.  To  overcome  this  difficulty,  we 

reverted  to  the  use  of  the  purely  repulsive  Mason-Vanderslice  exponential 

17  9 

Interaction  potential,  '  ' 

V  »  A  exp(-r/£).  (3.1) 

Then  the  dumb-bell  model  is  defined  as 


v_  -  V 

3  XM 

for  VXM 

>VYM 

and 

V  ■  V 

V3  YM 

for  Vm 

<VYM 

(3.2) 


This  change  should  not  affect  the  results  which  are  of  concern  here  since 

we  consider  only  dissociation  rate  measurements  at  high  temperatures 
o  3-5 

(2500  -  5000  K),  and  it  has  been  shown  that  the  potential  well  in  V, 
has  negligible  effects  on  reaction  rate  coefficients  at  such  elevated 
temperatures.  Table  I  summarizes  the  parameters  A^  and  used  in 
the  present  calculations. 


IIIB.  Reaction  Rates 

Table  II  shows  the  distribution  of  the  three  hundred  trajectories 
sampled  for  the  case  of  Ar  +  HC2.  at  2500  °K  (kT/D  ■  .0466)  with  e  ■  0, 
classified  according  to  the  nature  of  their  end  states  (see  Section  IIA). 
The  implications  of  Table  II  have  already  been  discussed  by  Keck.6  Of 
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course,  as  pointed  out  before,  his  empirical  correlation  formula  for 

(N/Nq)  does  not  apply  to  the  present  case.  In  addition,  Keck's  Table 

I  for  homonuclear  molecules  showed  entries  of  comparable  magnitude  for 

N(f|njf)  and  N(bjnjb);  whereas,  we  foind  that  N(jf|n|f)  »  N(b|njb)  for 

Ar  +  HCJt.  Hie  reason  for  this  difference  will  be  explained  in  Section 

IV  when  we  discuss  details  of  the  Ar  +  HC&  trajectories. 

The  recrossing  connection  factor  N/Nq  for  the  various  temperatures 

and  surfaces  sampled  are  plotted  versus  C  in  Figure  1.  No  definite 

temperature  dependence  can  be  determined  from  the  plot  (Keck's  previous 
2  6 

results  '  also  did  not  show  any  temperature  effect  on  N/Nq).  Although 

there  is  considerable  scatter  in  the  data,  the  least-squares  fit  line 

does  give  a  reasonable  representation  of  all  the  data  points.  The 

parameters  p  -  0.26,  b  ■  1  -  p  ■  0.74,  and  C  -  [N(0)/Nq(0)]  -  0.27  are 

all  obtained  from  this  fitted  line  (see  Equations  (2.19)  and  (2.21)). 

The  various  moments,  K  (c),  as  defined  by  Equation  (2.14)  were 

®1 

then  calculated  using  Equation  (2.15)  and  tabulated  in  Table  III.  The 
ratio  2K^/K2  «  a  (see  Equations  (2.26)  and  (2.27))  is  also  included  in 
the  table. 

In  order  to  complete  the  determination  of  the  transition  kernel, 

R(elt  ef ),  and  produce  a  map  of  its  relative  numerical  distribution,  we 

must  first  find  the  pre-exponential  function  G(e^)  (see  Section  IIB). 

2 

Keck  found  that  the  function 

G(e)  -  [1  -  (-e/6)1/2]/[l  +  (-e/6) 1/2 )  (3.3) 

was  able  to  represent  his  numerical  calculations  of&v(e)/^(0)  with 
sufficient  accuracy,  and  we  have  used  this  form  for  G(e^).  The  resulting 
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numerical  distribution  of  K^,  £f)/^(0)  is  show  in  Figure  2,  again  for 
Ar  +  HC l  at  2500  °K.  The  figures  in  the  boxes  represent  the  mean  value 
of  R(e^,  £f)/#^(0)  for  the  box.  Because  of  the  sampling  technique  employed, 
each  value  of  £  used  in  the  calculations  gives  data  in  a  quadrant  that 
touches  the  diagonal  -  e^,  at  a  value  of  (e^  +  e^)/2  ■  e  ■  e.  Data  in 

overlapping  quadrants  agreed  within  the  statistical  errors  and  were  averaged 
to  obtain  the  results  shown.  The  total  number  of  trajectories  sampled 
was  approximately  three  hundred  for  each  value  of  e(e  -  0,  -  1,  -2,  -3). 

The  results  shown  were  calculated  for  £^  >  (recombination);  however, 
since  R(eit  £f)  is  symnetrical  in  and  ef,  the  values  for  £^  <  ef,  may 
be  obtained  by  reflection  in  the  diagonal  ■  £,.  As  was  pointed  out  in 
Section  1IB,  the  map  shows  that  RCc^,  e^)  exhibits  a  sharp  maximum  ridge 
along  ■  ef,  which  corresponds  to  zero  energy  transfer;  this  dies  off 
rapidly  as  ^  increases.  Combining  Equations  (2.9)  and  (2.22),  we 

obtain 

R(eif  ef)^(0)  -  C(a2  -  b2/4)  G(ef)  exp(aef  -  Bei)  .  (3.4) 

Since  all  the  parameters  in  Equation  (3.2)  have  been  obtained,  we  can  use 
them  to  generate  the  function  R(£^,  £f)/<ft^(0).  The  parameter  a  used  is 
the  average  of  the  four  values  of  21^/^  in  Table  III.  The  result  for 
T  ■  2500  °K  is  plotted  in  the  upper  left  comer  of  Figure  2  as  a  contour 
map.  We  see  that  this  map  reproduces  all  the  important  features  of  the 
actual  distribution,  and  the  agreement  in  general  is  very  good. 

In  Figure  3  we  have  plotted  [R(A,  e)/^(0)  ]/[<R(e)/ft(0)  ]  versus  A. 

The  factor  [(R(e)/|ft(0) ]  is  included  for  normalization  purposes.  The  solid 
line  passing  through  the  bulk  of  the  data  points  has  a  slope,  a  ■  1.79, 


as  determined  above  (see  Table  Ill).  It  can  be  seen  that  for  A  <  0.5, 
the  data  points  deviate  considerably  from  this  line;  however,  the  region 
around  4*1,  has  the  most  important  contribution  in  the  determination 
of  the  nonequilibrium  correction  factor,  and  the  line  fits  the  data  very 
well  in  this  region. 

We  therefore  conclude  that  our  functional  form  for  R(e^,  e^)  and 
the  associated  parameters  determined  above  lead  to  a  very  good  representa¬ 
tion  of  the  transition  kernel.  At  the  same  time,  it  is  still  simple  enough 
so  that  the  steady-state  master  equation  can  be  solved  analytically  with¬ 
out  the  usual  reliance  on  a  great  deal  of  numerical  computation. 

The  nonequilibrium  correction  factors  may  now  be  calculated  from 
Equations  (2.13)  and  (2.36).  The  barrier  rate  coefficients  are  taken  from 
Reference  4  and  the  results  are  summarized  in  Table  IV  for  T  ■  5000°, 

3500°,  and  2500  °K.  It  should  be  noted  that  for  the  results  shown,  the 
corresponding  characteristic  energy  transfer  in  a  collision  is  between 
0.4  kT  and  0.6  kT.  The  statistical  errors,  based  on  one  standard  devia¬ 
tion,  are  estimated  to  be  about  10  per  cent  for  (N/Nq)  and  20  per  cent 

for  (k/k  ). 
e 

Figure  4  is  an  Arrhenius  plot  of  the  measured  (full  lines  and  open 

points)  and  calculated  (solid  points)  dissociation  rate  coefficients  of 

HCi  in  Ar.  It  shows  that  the  prasent  theoretical  predictions  are  in 

excellent  agreement  both  in  terms  of  magnitude  and  temperature  dependence 

with  the  experimental  measurements  of  Seery  and  Bowman,^  Fishburae,^ 

12 

and  Jacobs,  Cohen,  and  Giedt.  We  should  point  out  that  the  major  tem¬ 
perature  dependence  of  k^  is  determined  by  the  barrier  rate  coefficient  k 


o.  oj 
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The  correction  factors  (N/N  )  and  (k/k  )  lower  the  rates  and  Introduce 

o  e 

an  additional  relatively  weak  temperature  dependence  to  the  rate  coeffi¬ 
cients.  Error  bars  on  the  theoretical  points  show  estimated  statisti¬ 
cal  errors  of  -  20  per  cent. 

IV.  DISCUSSION 

13 

The  Landau-Teller  theory  of  vibrational  relaxation,  as  extended 

14 

by  Schwartz,  Slawsky,  and  Herzfeld  (SSH)  ,  which  considered  vibrational- 

translational  energy  transfer,  provided  both  a  successful  qualitative 

explanation  and  quantitative  correlation  of  the  vibrational  relaxation 

rates  of  most  of  the  homonuclear  diatomic  molecules  which  have  been 

investigated.  However,  when  applied  to  highly  asymmetric  heteronuclear 

diatomic  molecules  (e.g.,  HF  and  HC&),  the  theory  predicted  rates ^ 

several  orders  of  magnitude  slower  than  those  observed  experimentally.1^”2^ 

This  discrepancy  led  to  the  suspicion  that  rotational  motion  may  be  important 

in  vibrational  relaxation.  Much  work  has  been  done  in  this  respect,  and 

it  appears  that  by  properly  accounting  for  the  rotational  motion  of  the 

molecules,  the  very  high  vibrational  relaxation  rates  for  molecules  such 

as  HF  and  HC&,  as  well  as  the  effects  of  deuteration,  can  now  be  explained 

and  correlated.  Detailed  discussions  and  references  on  the  topic  of  rcta- 

21  2? 

tion-vib ration  coupling  can  be  found  in  papers  by  Moore,  Chen  and  Moore, 

23 

and  in  a  recent  comment  by  Nikitin. 

According  to  the  SSH  theory,  the  efficiency  of  energy  transfer  for 
a  collision  depends  on  an  "adiabaticity  factor,"  f(wT),  where  to  is  the 
vibrational  frequency  of  the  molecule  and  T  is  the  "effective  time"  of 
the  collision,  i.e.,  the  time  available  for  energy  transfer.  In  the 
theory,  the  relative  velocity  between  the  center  of  mass  of  the  molecule 
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1/2 

and  the  collision  partner  is  used  to  estimate  T;  thus  TOtO^  M)  * 

Wher*  UXY,H  ■  "XV 

collision  complex.  Such  reasoning  as  this  led  Keck  to  seek  a  correla¬ 
tion  of  the  recrossing  correction  factors  (N/N  )  for  homonudear  molecules 

o 

and  various  collision  partners  on  the  basis  of  th?  adiabaticity  parameter, 
1/2 

w  T  a  *  However»  *or  such  systems  as  HCJt  +  Ar,  the  hydro¬ 

gen  atom  moves  much  more  rapidly  than  does  either  the  C l  or  Ar  atoms  due 
to  its  much  lighter  mass.  As  can  be  seen  in  the  trajectories  shown  in 
Figure  5,  the  essential  motion  of  the  H  atom  is  that  of  rotation  about 
the  C i.  atom.  Thus  a  more  correct  collision  time  for  these  heteronudear 
systems  is  determined  by  the  motion  of  the  H  atom  relative  to  the  Ar  atom 
and  not  by  the  relative  motion  of  the  HCi  and  Ar.  Consequently,  the 
quantity  yHcJ^  r  in  the  SSH  theory  should  more  properly  be  replaced  by 

**HAr  “  “h  TBAr^mH  +  mAr^  88  "h*  t*ie  re^uce<*  10838  of  the  H-Ar  collision. 

This  point  is  clearly  illustrated  by  the  three  reactive  trajectories  we 
plotted  in  Figure  5.  The  duration  of  the  strong  interaction,  where  most 
of  the  energy  and  momentum  transfers  occur,  is  indicated  by  the  thicker 
lines  on  the  trajectories.  We  see  that  during  this  period  of  strong 
interaction  (corresponding  to  the  effective  collision  time),  the  Ar  atom 
hardly  moves  relative  to  the  C&  atom  whereas,  by  rotation,  the  H  atom 
passes  completely  through  the  reaction  zone.  Since  the  adiabaticity 

g 

factor  f(fa)T)  has  a  dependence  given  by  the  factor  exp(-  a  fa)  t)  ,  where 
both  a  and  3  are  of  order  unity,  the  change  from  l  „  to  increases 

AY  jfl  AN 

the  magnitude  of  the  adiabaticity  parameter  to  order  unity  for  such  sys¬ 
tems  as  HC&  +  Ar  and  thus  raises  the  calculated  vibrational  relaxation 


(m^  +  b^)  is  the  effective  reduced  mass  of  the 
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races  by  several  orders  of  magnitude  In  agreement  with  experimental 
observations.  Additional  trajectory  calculations  for  the  hetemuclear 
molecules  using  different  collision  partners  are  required  before  we  are 
in  a  position  to  deduce  correlation  formulas  for  the  correction  factors 
similar  to  those  obtained  by  Keck. 

The  inner  dashed  circles  in  Figure  5  indicate  the  equilibrium 
separation  for  HC&  molecules,  the  outer  dashed  circles  are  the  posi¬ 
tions,  z^,  of  the  rotational  maximum  evaluated  at  the  end  of  the  colli¬ 
sion.  In  general  the  following  observations  can  be  made  about  the  HCf.-Ar 
collisions: 

1.  For  non-reactive  collisions  (unlike  those  shown  in  Figure  5)  the  Ar 
atom  acts  as  a  block  which  prevents  the  H  atom  from  reaching  the  Cl 
atom  by  virtue  of  its  position  between  the  Cl  and  H  atoms;  i.e., 
rHCS,  >  rHAr*  These  collisions  are,  in  fact,  automatically  excluded 
from  our  samples  since  they  never  reach  the  reaction  zone. 

2.  On  the  other  hand,  participation  of  Ar  in  the  collision  is  necessary 

to  take  away  the  excess  energy  from  H*C&  aid  thereby  effect  a  recombina¬ 
tion.  Thus  typical  reactive  trajectories,  such  as  those  shown  in 
Figure  5,  may  be  divided  into  three  different  phases: 

(i)  The  H  atom  approaches  the  C&  atom  from  the  opposite  side  to 
the  Ar  atom  with  little  or  no  interference  from  Ar. 

(ii)  Due  to  the  action  of  the  strong  attractive  valence  forces 

between  H  and  Cl,  the  H  atom  cuts  into  a  line  joining  Ar  and 
Cl  and  transfers  energy  and  momentum  to  Ar.  This  is  the  zone 
of  strong  interaction  (thicker  line  on  the  trajectories)  during 
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whlch  most  of  the  energy  and  momentum  exchanges  are  effected. 
We  note  that  at  the  beginning  of  this  reaction  zone,  the  Ar 
atom  is  usually  very  near  the  turning  point  with  respect  to 
the  Cl  atom  so  that  its  kinetic  energy  is  very  low, 

(ill)  The  H  atom,  having  lost  sufficient  energy,  becomes  bound  to 
Ci.  atom  and  then  proceeds  to  vibrate  and  rotate  about  the 
latter;  whereas,  the  Ar  atom  moves  away  from  HC£  molecule 
with  an  increased  velocity. 

We  should  point  out  that  for  convenience  of  representation  in  two 
dimensions,  the  trajectories  shown  in  Figure  5  have  all  three  particles 
moving  in  the  same  plane.  However,  our  calculations  were,  in  general, 
three  dimensional. 

In  Section  IIIB  we  noted  that  N(f|n|f)  »  N(b|n|b)  for  Ar  +  HCf, 
collisions;  whereas,  Keck6  obtained  N(f|n|f)«  N(b|n|b)  for  all  the 
homonuclear  molecules  he  investigated.  The  reason  for  this  difference 
can  now  be  explained.  The  variational  reaction  rate  calculations  show 
that  the  contribution  to  the  reaction  rate  from  collisions  with  the  Cl 
end  of  the  molecule  is  negligible.  In  addition,  both  Ar  and  Cl  can  be 
regarded  as  almost  stationary  while  the  hydrogen  atom  moves  about  at 
high  average  velocities.  Therefore,  most  of  the  energy  exchanges  that 
may  change  the  state  of  HC£  must  come  from  the  li-Ar  collision,  with  very 
little  direct  influence  from  the  Cl  atom.  Our  sampling  procedure  ensures 
that  the  molecule  must  be  on  the  energy  surface  H  -  B  *  0  at  a  certain 
time  t  =  tQ,  and  the  trajectory  must  cross  that  surface  at  tQ  with  a 
negative  slope  (for  recombination);  i.e. ,  the  molecule  must  be  in  a  free 
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scate  (H  -  B  >  0)  at  some  time  t  <  tQ  and  must  be  losing  energy  to  the  Ar 
atom.  Consider  the  case  in  which  the  H  atom  was  initially  bound  to  the 
C&  atom  at  t  «  tfl.  The  result  of  the  H  +  Ar  collision  would  be  either 
(a)  a  net  loss  of  energy  from  H  to  Ar,  with  the  HC&  molecule  remaining 
bound,  or  (b)  a  net  gain  of  energy  by  H  from  Ar,  making  H  free  from  C£. 
Obviously  trajectories  (a)  are  not  sampled  by  our  procedure  since  they 
do  not  cross  the  surface  H  -  B  ■  0;  nor  do  trajectories  (b)  cross  the 
surface  in  the  desired  direction  since,  once  the  H  atom  becomes  free 
after  a  collision  with  the  Ar  atom,  its  rapid  motion  will  carry  it  away 
from  the  C&  +  Ar  system.  Thus  no  further  significant  energy  transfer 
between  H  and  Ar,  which  would  change  the  state  of  HC&,  is  probable; 
i.e.,  HC&  remains  free, and  those  trajectories  are  not  sampled.  On  the 
other  hand,  it  is  easy  to  see  that  if  H  was  initially  free  at  t  «  tQ, 
its  collision  with  Ar  would  produce  f  b  trajectories  that  satisfy  our 
sampling  requirements.  Consequently,  the  result  N(f|n|f)  »  N(b|n|b) 
for  Ar  +  HCS.  collisions  is  expected. 

For  homonuclear  molecules,  the  outcome  of  case  (b),  above,  may  be 
quite  different.  The  b  ■*>  f  transition  would  be  a  result  of  the  collision 
between  one  end  of  the  molecule  and  the  third  body,  as  before.  However, 
rotation  of  the  molecule  may  bring  the  other  end  into  a  collision  with 
the  third  body  and  could  thereby  produce  a  f  b  transition.  These 
trajectories  (b  f  ■*  b)  satisfy  our  sampling  requirements,  and  N(f|n|f) 
and  N(b|n|b)  may  therefore  have  comparable  magnitudes.  Typical  trajectories 
of  this  kind  which  illustrates  our  point  will  be  presented  when  we  report 
our  work  on  homonuclear  molecules. 


-20- 


It  should  be  noted  that  the  features  of  the  collision  mechanics 
discussed  above  would  be  suppressed  If  only  collnear  collisions  were 
considered.  Thus  the  one-dimensional  approximation  sometimes  used  in 
theoretical  analyses  of  this  nature  presents  a  picture  which  is  very  differ¬ 
ent  from  what  actually  is  the  case  for  both  homonudear  and  hetemudear 
molecules. 

For  a  collision  to  be  reactive  (consider  recombination),  the 
molecule  must  start  in  a  free  state  and  end  in  a  bound  state.  The 
corresponding  energy  changes  must  be  such  that  H  -  B  <  0,  where  H  is 
the  Hamiltonian  of  the  molecule,  and  B  is  the  height  of  the  rotational 
barrier.  Figure  6  shows  a  distribution  plot  of  AH  andAB  for  101  reactive 
trajectories  out  of  the  three  hundred  trajectories  sampled  for  HCl  +  Ar 
at  2SC0  °K  on  the  surface  S(0).  It  is  apparent  that  recombination  (the 
change  from  H-B>0toH-B<0)  is  achieved  by  two  equally  Important 
mechanisms,  i.e.,  reducing  H  and/or  increasing  B.  Our  work  on  homonudear 
molecules  (to  be  reported  separately),  on  the  other  hand,  shows  that  most 
recombinations  are  achieved  by  reducing  H.  This  furnishes  further  evidence 
that  rotational  energy  transfer  in  asymmetric  molecules  pl  ays  a  relatively 
important  part  in  chemical  reactions.  The  energy  transfers  AH  and  AB  for 
the  three  trajectories  shown  ir.  Figure  5  are  indicated  by  the  filled  points 
in  Figure  6.  They  represent  different  ratios  of  AB/AH. 

The  calculations  presented  in  this  paper  are  for  one  system  (Ar  +  HCS.) 
and  over  a  limited  (high)  temperature  range  (2500°  -  5000  °K) .  Work  on 
other  diatomic  molecules  is  in  progress  which  will  cover  the  high  tempera¬ 
ture  dissociation  region  and  the  low  temperature  recombination  region.  A 
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number  of  Interesting  questions  can  be  answered  by  such  calculations; 
e.g.,  what  Is  the  physical  Interpretation  of  the  parameters  a  and  8 
Involved  In  the  equilibrium  transition  rates,  and  what  will  be  the  con¬ 
tribution  from  the  "complex"  mechanism  to  the  overall  recombination 
rate  coefficient  at  low  temperatures?  In  any  event,  the  agreement  between 
the  theoretically  predicted  dissociation  rate  coefficients  and  experimental 
measurements  for  the  case  of  HCS,  +  Ar  presented  here  undoubtedly  provides 
another  successful  test  for  the  classical  modified  phase-space  theory 


of  reaction  rates. 
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TABLE  I.  Maaon-Vanderslice  Potential  Parameters 


o 


System 

A(eV) 

KA) 

Ar  +  H 

268 

0.35 

Ar  +  Cl 

1060 

0.25 

iffimMMk  wk- 
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TABLE  II.  Distribution  of  trajectories  with  respect  to  class  of 

reaction  and  number  of  traversals  of  the  barrier  surface 
(e  -  0),  for  Ar  +  HCS.  at  2500  °K(kT/D  «  0.047) 


Class 

n 

N 

Total 

1 

2 

3 

(f |n|b) 

91 

8 

99 

(f|n|f) 

142 

9 

2 

153 

(b  |  n  |  b ) 

10 

4 

14 

(b  |  n  |  f ) 

16 

2 

18 

Total 

284 

Recrossing  correction  factor 

(N/N  )  -(91  +  8/2)/284  -  0.33 
o 
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TABLE  III.  Momenta  [K  (C)]  for  Ar  +  HC& 

III 


T(°K) 

kT/D 

c 

*0 

B 

K2 

2Kx/K2 

5000 

.093 

0 

1.09 

.245 

.159 

3.08 

-1 

.501 

.182 

.158 

2.30 

-2 

.505 

.102 

.079 

2.58 

-3 

.348 

.080 

.043 

3.72 

mean 

2.92 

3500 

.065 

0 

1.34 

.233 

.174 

2.68 

-1 

.693 

.190 

.242 

1.57 

-2 

.321 

.130 

.130 

2.00 

-3 

.464 

.112 

.095 

2.36 

mean 

2.15 

2500 

.0A7 

0 

1.71 

.331 

.332 

1.99 

-1 

.675 

.192 

.247 

1.56 

-2 

.555 

.202 

.241 

1.68 

-3 

.364 

.157 

.162 

1.94 

mean 


1.79 
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TABLE  IV.  Dissociation  Rate  Coefficients  for  Ar  +  HCfl, 


T(°K) 

Hi 

-*|Q 

a 

N 

N 

o 

b 

a 

k 

k 

e 

kd 

5000 

.093 

5.69”13 

.27 

.74 

2.92 

.30 

4.55‘U 

3500 

.065 

7.26"15 

.27 

.74 

2.15 

.42 

8.27"16 

2500 

.047 

1.81"17 

.27 

.74 

1.79 

.50 

2.44-18 

k®  and 
a 

k ,  have 
a 

units  cm3 

molecule 

1 

sec 

1 

-28- 


Figure  1 

Figure  2 

Figure  3 

Figure  4 

Figure  5 
Figure  6 


FIGURE  CAPTIONS 


Recrossing  correction  factors,  N/N  versus  energy  of  the 
surface  from  which  samples  were  tafcen.  #  ,  T  -  5000  °K; 

0 ,  T  -  3500  °K;  O .  T  -  2500  °K.  Solid  line  is  the  best 
least-squares  fit  for  all  points.  The  slope  of  this  line 
determines  the  parameter  b. 

Map  of  numerical  distribution  of  the  function  R(£j.,  £f)4Rv(0), 
where  R(£j.,  £f)  is  the  equilibrium  transition  rate  andftv(0) 
is  the  variational  rate  for  the  surface  e  ■  0.  The  contour 
map,  calculated  from  Equation  (3.2)  using  the  fitted  parameters, 
is  plotted  in  the  upper  left  corner  to  avoid  overcrowding  the 
figure.  Comparison  may  be  done  by  reflection  through  the 
diagonal  £^  ■  £^. 

The  function  [R(A,  e)/0?v(O)  ]/ El8(e)/lft(0)  ]  is  plotted  versus  A, 
the  absolute  value  of  energy  transfer.  £  ■  0;  O*  e  ■  -1; 

A,  £  ■  -2;  e  ■  -3.  The  factor  [02(£)/A(O)]  is  included 
for  normalization  purposes  so  that  the  data  can  be  better 
represented  by  a  single  line.  The  solid  line  which  passes 
through  the  bulk  of  the  data  has  a  slope  a  taken  from  Table  III. 
See  text,  Section  IIIB. 

Comparison  between  theoretical  predictions  of  the  dissociation 
rate  coefficient  and  experimental  measurements,  as  a  function 
of  T.  ♦  ,  theoretical  predictions;  O.  Seery  et  al;  (F), 
Flshbume;  (J),  Jacobs  et  al;  dashed  curve  represents  the 
theoretical  barrier  rate  coefficients. 

Typical  reactive  trajectories.  The  thicker  lines  on  the 
trajectories  indicate  zones  of  strong  interaction. 

Distribution  plot  of  reactive  trajectories  in  the  [AH,  AB] 
plane. 
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DISSOCIATION  RATE  COEFFICIENT,  kd  cm3  molecule'1  sec' 
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